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the impurities makes the fluctuating amplitude mode visible in the optical response. 
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1 Introduction 



Quasi one-dimensional conductors like the transition metal chalcogenides are charac- 
terized by a nested Fermi surface. This renders them unstable to the Peierls transition 
when electron-phonon backscattering between the two sheets of the Fermi surface 
become relevant. We call this a Peierls-system (PS). The ensuing charge density 
wave (CDW) which develops below a transition temperature Tp of the order of 200K 
shows many unusual properties, especially in its electronic responses [[j]. 
In recent years it became increasingly evident that the electronic properties of PS 
in their normal phase (T > Tp) are also unusual. Photoemission studies (0, H and 
earlier work cited therein) point towards possible non-Fermi-liquid (NFL) behaviour. 
The microwave and optical response also deviates from the Drude predictions for a 
normal metal j|, Specifically, the real part of the complex conductivity function 
-Recx(u;, T ) for the chain direction and at room temperture T shows the presence 
of a pseudo-gap at u « 2A where A is the zero temperature half gap of the 
corresponding CDW. In addition, a peak structure appears in the far infrared well 
below the pseudo-gap. This peak resembles the pinned Frohlich mode at about the 
same frequency but seen in the fully developed CDW. Similar results were also found 
in PS films @. 

In principle these features are known for a long time from the CDW material K 2 Pt{C TV^-Bro^- 
3(H20) (KCP) in which fluctuation effects due to intrinsic disorder are prominent. 
In |7| the ac-conductivity of KCP is modelled by a dielectric function which takes 
pinning-unpinning fluctuations into account. 

A related explanation for the optical data is provided by the concept of fluctuating 
CDW segments which behave similar to the fully developed CDW. These segments 
are pinned by random impurities which break the phase invariance of the equation 
for the fluctuating order parameter. 

In view of the one-dimensional nature of fluctuations above the temperature T* > 
Tp below which transverse fluctuations set in to initiate the phase transition, the 
observed effects could also be a signature of NFL behaviour. The latter seems to be 
established for some of the Bechgaard salts |], |^, ID] where CDW fluctuations are 
not important. 

In the case of PS, the attractive backscattering and the softening of the Peierls 
phonons make a plain Luttinger-liquid scenario unlikely. Voit |11], advocates a 
Luther-Emery state [13| for the electrons in blue bronze as explanation of the pho- 
toemission spectra. In this case a spin gap would open before true CDW formation. 
Such conclusions, however, are not undisputed. Shannon and Joyce |14| argue in 



favour of a model [T5|, [16j for a fluctuating Peierls system (FPS). The unusual plas- 
mon dispersion in quasi one-dimensional conductors can also be understood in terms 
of a conventional band picture flI7|. There seems to be good reason to persue the 



FPS concept for one-dimensional metals with CDW fluctuations. 

The FPS model has recently been extended by McKenzie [IS]. He calculates the 



one-particle Green's function for Gaussian order parameter fluctuations and large 
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correlation length using Sadovskii's exact method JT9|. He points out that the elec- 
tron spectral function of FPS is of NFL type as found earlier in the same context [ 20 1 . 
McKenzie also calculates renormalized coefficients for Ginzburg-Landau functionals 
of PS. The model which we will solve for the collective conductivity corresponds to 
the FPS concept. 

The present paper computes the frequency dependent conductivity from one-dimen- 
sional conductive order parameter fluctuations using a modified, linear Ginzburg- 
Landau-Langevin (GLL) equation for CDW including phase breaking by impurities. 
In superconductivity (SC) where phase pinning does not exist, this part of the con- 
ductivity was first investigated by Aslamazov and Larkin [21| (AL). The CDW version 
of the AL-theory was given in |22| . 



From SC it is known |23|, [24]] that there are two further contributions from order 



parameter fluctuations to the dc conductivity: A resistive contribution from the 
reduction of the single particle density of states which is related to the pseudo gap 



and the anomalous Maki-Thompson p5[, ^g] term which is conductive. In CDW the 



latter becomes resistive J^J and dominates the dc conductivity near the transition. 
The issue of one-dimensional collective dc conductivity in CDW was strongly de- 



bated in the seventies |28], ^y, [22, 30, [31, p2| . In |28| the idea of paraconductivity 



in one-dimensional metals with dominant electron-phonon interaction was advanced 
using a GL approach. This result was criticized in f29| where pinned collective fluc- 
tuations were shown to reduce the dc conductivity. These authors also studied the 
corresponding results for the Hubbard and the Tomonaga-Luttinger model. Detailed 
microscopic studies of FPS in P2L |5(J find both resistive and conductive fluctuations 
but neglect phase pinning. The dc paraconductivity of commensurate FPS was stud- 
ied in [^| in a GL context. This paper served as a starting point for the present 
work. 

The paper is organized as follows: Sec. 2 develops the GLL approach for CDW and 
derives the known results for unpinned conductive fluctuations. Sec. 3 introduces a 
modified GLL-equation where phase invariance is broken in a way which simulates 
local pinning by random impurities, and evaluates the basic correlation function of 
the order parameter fluctuations. This result is used to calculate the frequency de- 
pendent collective conductivity exactly within the model. The complicated formula 
is evaluated approximately for pinning frequencies small in comparison to the fre- 
quency of amplitude fluctuations in Sec. 4. Two appendices present mathematical 
details. 



2 Ginzburg-Landau— Langevin Approach 



It is interesting to formulate the problem of fluctuation ac-conductivity in CDW 
without pinning using the GLL-method. One expects to find close similarities to 
the elegant formulations for SC |33|, [34[]. However, it turns out that one must go 



beyond overdamped dynamics which in case of CDW would only give an instanta- 
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neous response in the current correlation function and thus a frequency independent 
conductivity. 

Our starting point is the GLL equation directly in terms of the gap fluctuations 



A fc (t) + 70 A fc (t) + u 2 k A fe (t) = r fc (t). 



The parameters 70 and uik are taken from 
spond to the rigid lattice values in 



(1) 

Static parameters below corre- 
For simplicity the renormalization of the 
rigid lattice (RL) values due to fluctuations as proposed in |1| is not considered. To- 
gether with Hartree-Fock corrections to the linear GL-equation as in [p5 |, it would 
extend the region of applicability of the linearized approach which is marginal at 
the RL-level. This and the absence of resistive fluctuations prohibit a quantitative 
comparison with experiments. Below the CDW transition temperature McKenzie's 



approach |18| is neccesary to give the characteristic optical absorption which has 



been measured in [3£ 



The damping constant 70 is 



where 



7o 



fm 
l 8k B T 



(2) 



(3) 



is the frequency of the amplitude mode of the fully developed CDW [j37| , ujq is the 
bare frequency of the 2£;^-phonon which goes soft, and A is the electron-phonon 
coupling constant. 

The actual frequency oj^ of the amplitude fluctuations in ([l]) is 



(1 + ee 



(4) 



with 



uj = uj A e RL . (5) 

For defmiteness we assume underdamping (70 < 2uq) which is the case in the linear 

fluctuation regime sufficiently above the transition temperature. 

Without a form of selection rule breaking neither uja nor ujq can be observed in 

optical conductivity measurements. The amplitude mode uja can be seen, however, 

in Raman scattering []38], |39|. The fluctuating amplitude mode above the transition 

temperature is observed in neutron scattering studies |40], 

The correlation length £ in (PS) is 



£ 2 — Co/ € RL, 



(6) 



with €rl given by 
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T 

e RL = \n—, (7) 

J-RL 

where Tm is the rigid lattice mean field transition temperature. The reference length 
£ is given by jH 



2 7Cm 2 v 2 F 

^° 167T 2 {k B T) 2 ' 

The complex Gaussian Langevin force T has zero mean. Its correlation function 



(T k (t)Tl,(0)} = (\A k \ 2 } 2 lo u; 2 k 5 k , k ,8(t) = 2 7o ^ 2 -§- 6 k , k , 5(t) = A8 kjk ,8(t). (9) 

Jo 

is constructed in such a way as to give the fluctuation intensity 

<|A fe | 2 > = ^. (10) 
a k 

The latter follows from the linear free energy functional 

F = J2a k \A k \ 2 , (11) 

k 

with 

a k = f (l + k 2 e), fo = ^- (12) 
Note that (|) impl ies a spatial correlation of the order parameter according to 

(AM)A*(0,0)> = ^/ ^e^(|A fc | 2 ) = (^^) e'^^^e^. (13) 

In the next step the one-dimensional conductivity is computed from the classical 
Kubo formula 

T. roc 

a (u) = -tL. / dte lujt (J(t)J(0)), (14) 



k B T Jo 

where L is the sample length. 

The collective current density was calculated in |32| [] and reads: 

j(x, t)=i h - (A(x, t) A*(z, t) - A(x, t) A*(x, t)). (15) 



1 S.N. Artemenko informed the author that he obtained the same expression for the collective 
current density by the Keldysh approach to CDW dynamics. 
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The collective current is proportional to the time derivative (f of the order parameter 
phase as for the fully developed CDW but the prefactor is different. The coefficient 



6 in (JT5D is ]3T 



e 



M B ThvbJ ' ^ 

involving the backward scattering rate i/f, due to random static scattering centers. 
This formula holds in the pure limit when the electron scattering rate obeys Hu = 
H"f + Vb/2) < 2nk B T. 

The homogeneous current density J in (fTEf) is related to j by 



;i7) 



i r L 

J(t) = — dxj(x,t) = j k= o(t). 
L Jo 



In the linear setting of (|l|) not only T obeys Gaussian statistics but also A and exact 
Gaussian decoupling gives 



(mm) = j e t) 2 - c (k, t)c(k, t) 



provided the correlation function 

C(fc,t) = (A fe (t)A*(0)) 
is real and even in t. C(k,t) is evaluated from ([!]) and explicitly given by 



7o 



C(M) = (|A fe | 2 ) exp(-^|t|) 



cos Dj~t + sin DM 
2D h 1 1 



with 



(19) 



(20) 



D u 



2 7o 



'UJ 



(21) 



This leads to 



C 2 - CC = (\A k \ 2 )l ul exp(- 7o |t|). (22) 

Note that all oscillating terms in the correlation functions cancel out leaving a purely 
relaxational response. If one were to use the correlation function 

C(M) = (|Afc| 2 >o exp(- 7fc |t|), (23) 

for the overdamped version of (jl]) with 7^ = ui^/jo one would get an instantaneous 
response (J(t) J(0)) oc 8(t) and hence a frequency independent conductivity. 
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Calculation of the conductivity using ( p2|) gives, however, the correct result given in 
3| 



Rea(u)= a F ^L^, (24) 
7o +^ 

irrespective of the relation between 70 and uo. The scale-value op of the fluctuation 
conductivity is 

L 2 A% 2 

° F ~ i«; t 2c 3' ( 25 ) 



and coincides with the result [RTJ. Explicitly o> reads 



a F = . (26) 

'7C(3) eR i 



— 1/2 

The conductivity shows the mean-field critical behaviour op oc e^/ . In the picture 
of a metal with order parameter fluctuations the collective conductivity adds to the 
normal state conductivity = Sefy f / (AnfiVb) . 

Formal calculations in higher spatial dimensions require a momentum cut-off in con- 
trast to SC. This is related to the form of the collective current density QTST). 

3 Breaking of Phase Invariance 

The space-time version of (0) is 

A(x, t) + 70 A(x, t) + oo 2 (\ - e A(x, t) = r(ar, t). (27) 

The simplest way to break the phase invariance of this equation is to add a pinning 
term 

2u 2 \A(x,t)\ cos<p(x,t), (28) 
to the left hand side which is a simple local coupling. This is clearly not the general 



starting point to treat pinning by random impurities f|3|. However, it will become 
evident later that this approach simulates local pinning because the final pinning 
frequency is proportional to the impurity concentration. 
To arrive at (|28|) we start from the more general form 



uj 2 s ]T h(x - Xl ) (A(xi, t) + A*( Xi , t)). (29) 



The impurities are locally coupled to the order parameter. The real structure function 
h(x) transmits the effect to the order parameter at x. The terms A*(xi, t) break phase 



7 



invariance by modelling backward scattering. The scale frequencies u s and Ui are 
different from the final pinning frequency. 

We make two further asumptions: The function h is a contact interaction h(x) = 
li5(x) with a scattering length l t . The crudest assumption is, however, 

J2^ n i / dx h (30) 



where n, is the density of impurities. This requires > 1 and amounts to an early 
impurity average. Introducing the scale frequency 

uf = u 2 s nili (31) 

then leads to (p8|) . 

This admittedly crude model has the advantage to allow for an exact solution. 
The complete GLL-equation which replaces (|l|) can be written: 

A k (t) + 7o A k (t) + ul A k (t) + ul ( A fc (t) + A*_ fe (t) ) = Y k (t) . (32) 

The complex order parameter A(x,t) is decomposed into real and imaginary parts 
U(x,t) = ReA(x,t) and V(x, t) = ImA(x,t) giving 

A k = U k + iV k , (33) 
with complex U k and V k which satisfy the usual reality conditions 

U- h = Ut, V- h = V;. (34) 
The Langevin equation splits into two equations for U k and V k : 

U k (t)+ l0 U k (t)+u 2 k U k (t) + 2^U k (t) = F Uk (t), (35) 
V k (t)+ l0 V k (t)+uj 2 k V k (t) = T Vk (t). 

Here Tu k (t) and Tvk(t) are the Fourier transforms of the real and imaginary part of 
T(x,t), respectively. It is possible that the impurities modify the thermal random 
force T(x,t). In our model we assume that this is not the case. This assumption 
is reasonable for uj,i uj q . The random forces T Uk and T Vk are then independent 
Langevin forces with the same statistical properties as T k (t) (cf. (|9|)) but only half 
its strength. The two equations (|35D become independent and are both isomorphic 
with (|IJ). However, the frequencies for the £/-modes are modified and change their 
fluctuation intensities: 



l^| 2 ) = (l + 24r 1 ^ (W k \ 2 ) = ^. (36) 
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Hence the intensity of the fluctuating order parameter is reduced: 



;o l + 2(ui/u k y 



lAfcl 2 ) = (|A fc | 2 ) ( 



(37) 



A thermodynamic derivation for (|36|) is given in Appendix A. In view of the real- 
istic condition Ui <C Uk the renormalization of the mean square order parameter is 
irrelevant. 

The order parameter correlation function becomes 



with weights 



and the replacement 



C(k,t) = P+ C+(k,t) 



P+ 2(l + 2(^))' P ~ 



1 

2' 



(3t 



(39) 



(+) 



><4 + M - 1J T 



(40) 



in the expression (BO) for C(k,t) in order to get C + (k,t) while C_(k,t) remains 
unchanged, i.e. formally = Z^, ^ = cjfc. This solves completely the GLL- 
equation (|32|) . 



4 Discussion of Fluctuation Conductivity 

We define the wave-number dependent pinning frequency 

Wp(fc) = ^f. (41) 

From (|31|) it is seen that uj p (k) is proportional to the linear impurity concentration. 
Thus our model simulates local pinning. 

Using the condition Ui <C D k the result of Appendix B leads to the following expres- 
sion for the real part of the fluctuating conductivity 



Rea(uj) 



L b 2 



E(IA fe | 2 )^ fc 2 



7o 
7o + ^ 



(42) 



AL 
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+ 



7o 



+ 



7o 



(3-7o 2 M 2 )(7o 2 + 4^)-^ 

4D 2 (7 2 +w 2 + 4jD 



2^2 



Even before the fc-summation is performed three different contributions to the fluc- 
tuation conductivity can be discriminated: the relic of the AL-conductivity (AL) 
centered at zero frequency, a pinned mode (P) near the frequency o> p (0), and a weak 
structure (A) associated with the fluctuating amplitude modes to^. The latter results 
from selection rule breaking by the impurities. Thus traces of the fluctuating ampli- 
tude mode should be seen in the optical conductivity. A similar case regarding the 
pinned Frhlich mode is found in the fully developed CDW [44, 45|. 
The /c-summation is easily done for the AL-part and gives 



1 7 2 
Re a{uj) AL = - a F 

2 % + u 2 



This is exactly half the result fl24]) . The spectral weight W according to 

IT = / dujReo(uj) 



is 



W AL 



71 



A lenghty but exact calculation gives for the P-mode 



Rea{u)p = a F [Ji(u) + 2Re Jz{oS)\, 



with 



1 



7 2 + ^ 2 - Aufhl 



j f \ -"-2 

~ 2 7 ° (c 2 + [ 7o + 2uf/ l0 } 2 ) ( w 2 + [ 7o - 2ufh } 2 ) 



and 



J 2 M = 2^ 



7o 4^ 4 / 7o 2 - ( 7o + iuf ^(4^2 _ 7q2) ( 7q ~ luJ f ~ A(jJ 4 



(43) 



(44) 



(45) 



(46) 



(47) 



(4* 



Note the non-algebraic structure of the conductivity due to J 2 . The spectral weight 
associated with Re <t(uj) p is independent of pinning parameters and given by 



7T 



W P = - 7o 0>. 



(49) 
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It adds the missing half to the total spectral weight Tcyoap of the unpinned fluctuation 
conductivity. 

The amplitude mode is treated approximately. Assuming 70 <C 2u and retaining 
the fc-dependence only in the prefactor it is found 



f To 2 (12ag ~ ^) 

64 F \ ( W 2 _ 4tu 2) 2 + 4 ^2 7 2 

The amplitude mode has a peak near 2u; - Its spectral weight is small and given 
by 

^ = (sS (51) 

Fig. 1 shows the pinned fluctuation conductivity in comparison with the unpinned 
case neglecting the weak amplitude mode. 



5 Summary 

The Ginzburg-Landau-Langevin method is developed for the fluctuation conductiv- 
ity in charge density wave systems above the transition temperature when fluctua- 
tions are one-dimensional. An additional phase breaking term due to impurities is 
introduced and its consequences for the fluctuation conductvity is evaluated. It is 
found that the spectral weight of the unpinned fluctuations is partly redistributed 
into a pinned mode around a pinning frequency in the far infrared, as seen in exper- 
iments. In addition, selection rule breaking by the impurities enables traces of the 
fluctuating amplitude mode to appear in the optical response. 
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A Appendix: Pinned Fluctuation Intensity 



The deterministic part of the GLL-equation ( |32] ) can be expressed as in terms of real 
variables x^ v as 



with the energy functional 



^ + 70 ^ = -2^^? (al) 



F = J2 a kA, (a2) 

k.v 
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with 



r a fc (l + 2&(*0), ^ = 1,2, 
akv ~\a k , z/ = 3,4, 



and 



^ 2 



&(*) = -\. (a4) 

We have used the decomposition ( |3~3"D and splitted the components U k and V k into 
real and imaginary parts: 

£4 = arjfei + ix k2 , 14 = x fc3 + ix k4 . (a5) 

Though in our model the random forces generate no dependence among the x ku , the 
latter are not independent since the x kv are even under k —>■ —k for v = 1, 3 and odd 
for v = 2,4. In terms of independent x kv the energy F becomes 

f = 2EE «***L = E E ( a6 ) 

fc>0 t/=l fc>0 v=\ 

The statistical average of the independent variables is simply 

/ 2 v J dx ku x 2 ku exp(-(3F ku ) 
{Xkv) fdx ku exp(-f3F ku ) ' 

and gives the results (|36|). 

B Appendix: Calculation of Fluctuation Conduc- 
tivity 

We use (p8D in the Kubo formula 

r A 2 poo r . -. 

a{uj) = k^f2Jo dte^J2[C(k } tf-C(k 7 t)C(k,t)]. (bl) 

Splitting the correlation function C(k,t) into its constituents gives 

C 2 - CC = ]T vl (pi - C V C V ) + P+ p_ {2C + C_ - C + C_ - C_CV} . (b2) 
The result (122[) translates into 
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Cl - C V C V = (\A k \ 2 ) 2 U U) ) 2 exp(- 7 o|t|), v = ±. 



(b3) 



Hence 



C 2 - CC = (\A k \ 2 ) 2 exp(- lo \t\)p 



+ oo 2 + - {2C+C- - C + C_ - C_(7 + } 



(b4) 



A somewhat tedious calculation gives a formally exact expression of the frequency 
dependent fluctuation conductivity: 



<t(u) 



L b 2 r°° 



k B T 2 Jo 



dtexp(iujt - 7 t) ^p+(|A fc | 2 )o 



(b5) 



,2 , 1 (\ u J i + 2u 2 u 2 - 1 2 {u 2 + u 2 )/A 



{cos{D^ - D^)t - cos{D^ + D^)t 
+((4 + u 2 ) (cos(D k +) - D^)t + cos(D[ +) + D[7 ] )t 



D 



(+) 



+ 



an(D^-DD\t\ 



(-)^ 



D 



(-) 



(+) 



(-) 



This result is too complicated to be discussed in full generality. We will take advan- 
tage of the fact that in practice the condition uof <C D 2 is fulfilled and perform an 
expansion of (b_5) with respect to 



D 2 ^ 



This gives 



(b6) 







B 



2(u 2 k + u 2 ) cos -^t + uf 



2Di 



cos 2D k t 



(b7) 



^sin^|t| + ^4sin2^|t| 
D k D k 2DI * M 



4 ^fe - 7o/ 2 



2^1 



cos — —t. 
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It is easier to do the time integral for the real part of the conductivity in the ap- 
proximated version of (p5D . The imaginary part follows from the Kramers-Kronig 
relation. The relevant terms up to order 0Jp{k) but neglecting small corrections of 
numerical constants of order uf are given as (|4^) in section 4. 
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Figure 1: Real part of scaled fluctuation conductivity (sum of eq. (43) and eq. 
(46)) in comparison to the unpinned case (gray line: eq. (24)) as function of scaled 
frequency. Conductivity unit is ap according to eq. (25) and frequency unit is the 
damping constant 7 (cf. eq. (2)) of the fluctuating amplitude mode. In this unit 
the following values were chosen: amplitude mode frequency ojq = 40 and pinning 
scale frequency Ui = 16. The actual pinning frequency is seen to be near 670. 
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